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Chiral Extrapolation of Hadronic Observables 
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One of the great challenges of lattice QCD is to produce unambiguous predictions for the properties of physical 
hadrons. We review recent progress with respect to a major barrier to achieving this goal, namely the fact that 
computation time currently limits us to large quark mass. Using insights from the study of the lattice data 
itself, together with the general constraints of chiral symmetry, we demonstrate that it is possible to extrapolate 
accurately and in an essentially model independent manner from the mass region where calculations will be 
performed within the next five years to the chiral limit. 
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1. INTRODUCTION 

Involving as it does a finite grid of space-time 
points, lattice QCD |Q,p| requires numerous ex- 
trapolations before one can compare with any 
measured hadron property. The continuum limit, 
a — > (with a the lattice spacing) , is typically un- 
der good control M. With improved quark and 
gluon actions the 0(a) errors can be eliminated 
so that the finite-a errors are quite small, even at 
a modest lattice spacing - say 0.1 fm (4|^|. 

In contrast, the infinite volume limit is much 
more difficult to implement as the volume, and 
hence the calculation time, scales like A 4 . This 
limit is also inextricably linked to the third ex- 
trapolation, namely the continuation to small 
quark masses (the "chiral extrapolation"). The 
reason is, of course, that chiral symmetry is spon- 
taneously broken in QCD, with the pion being a 
massless Goldstone boson in the chiral limit. As 
the lattice volume must contain the pion cloud of 
whatever hadron is under study, one expects that 
the box size, L, should be at least 4m" 1 . At the 
physical pion mass this is a box 5.6 fm on a side, 
or a 56 4 lattice with a = 0.1 fm. This is roughly 
2 4 times as big as the lattices currently in use. 
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Because the time for calculations with dynami- 
cal fermions (i.e. including quark-anti-quark cre- 
ation and annihilation in the vacuum) scales as 
m~ 3 ' 6 ||, current calculations have been limited 
to light quark masses 6-10 times larger than the 
physical ones. With the next generation of su- 
percomputers, around 10 Teraflops, it should be 
possible to get as low as 2-3 times the physical 
quark mass, but to actually reach that goal on an 
acceptable volume will require at least 500 Ter- 
aflops. This is 10-20 years away. 

Since a major motivation for lattice QCD must 
be to unambiguously compare the calculations of 
hadron properties with experiment, this is some- 
what disappointing. The only remedy for the 
next decade at least is to find a way to extrap- 
olate masses, form-factors, and so on, calculated 
at a range of masses considerably larger than the 
physical ones, to the chiral limit. In an effort to 
avoid theoretical bias this has usually been done 
through low-order polynomial fits as a function 
of quark mass. Unfortunately, as we discuss in 
sect. 2, this is incorrect and can yield quite mis- 
leading results because of the Goldstone nature 
of the pion. 

Once chiral symmetry is spontaneously broken, 
as we have known for decades that it must be in 
QCD and as confirmed in lattice calculations, all 
hadron properties receive contributions involving 
Goldstone boson loops. These loops inevitably 
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lead to results that depend on either logarithms 
or odd powers of the pion mass. The Gell-Mann- 
Oakes-Renner relation, however, implies that 
is proportional to the square root of m q , so loga- 
rithms and odd powers of are non-analytic 
in the quark mass 0], with a branch point at 
TTtq = 0. One simply cannot make a power se- 
ries expansion about a branch point. 

On totally general grounds, one is therefore 
compelled to incorporate the non-analyticity into 
any extrapolation procedure. The classical ap- 
proach to this problem is chiral perturbation the- 
ory, an effective field built upon the symmetries 
of QCD ||. There is considerable evidence that 
the scale naturally associated with chiral symme- 
try breaking in QCD, A x sb, is of order 47r/ T , or 
about 1 GeV. Chiral perturbation theory (xpt) 
then leads to an expansion in powers of m n /A x sB 
and p/A x sb, with p a typical momentum scale 
for the process under consideration. At 0(p 4 ), 
the corresponding effective Lagrangian has only a 
small number of unknown coefficients which can 
be determined from experiment. On the other 
hand, at Q(p 6 ) there are more than 100 unknown 
parameters |10| , far too many to determine phe- 
nomenologically. 

1.1. Convergence of xpt ? 

Another complication, not often discussed, is 
that there is yet another mass scale entering the 
study of nucleon (and other baryon) structure 
|TT| This scale is the inverse of the size of 
the nucleon, A ~ Since A is naturally more 

like a few hundred MeV, rather than a GeV, the 
natural expansion parameter, m^/A, is of order 
unity for ~ 2 — 3mP hys - the lowest mass 
scale at which lattice data exists. This is much 
larger than m„-/A x sB ~ 0.3 — 0.4, which might 
have given one some hope for convergence. As it 
is, the large values of m^/A at which lattice data 
exist make any chance of a reliable expansion in 
traditional (dimensionally-regulated) chiral per- 
turbation theory (xpt) fairly minimal fll4|| . Even 
though one has reason to doubt the practical util- 
ity of xpt, the lattice data itself does give us some 
valuable hints as to how the dilemma might be re- 
solved. The key is to realize that, even though the 
masses may be large, one is actually studying the 



properties of QCD, not a model. In particular, 
one can use the behaviour of hadron properties 
as a function of mass to obtain valuable new in- 
sights into hadron structure. 

1.2. Where the Constituent Quark Picture 
Might Work 

The first thing that stands out, once one views 
the data as a whole, is just how smoothly every 
hadron property behaves in the region of large 
quark mass. In fact, baryon masses behave like 
a + bm q , magnetic moments like (c + diriq)" 1 , 
charge radii squared like (e + /m g ) _1 and so on. 
Thus, if one defined a light "constituent quark 
mass" as M = M + cm q (with c ~ 1), one would 
find baryon masses proportional to M (times the 
number of u and d quarks), magnetic moments 
proportional to M~ x and so on - just as in the 
constituent quark picture. There is little or no 
evidence for the rapid, non-linearity associated 
with the branch cuts created by Goldstone boson 
loops. Indeed, there is little evidence for a sta- 
tistically significant difference between properties 
calculated in quenched versus full QCD! How can 
this be? 

The natural answer is readily found in the ad- 
ditional scale, A ~ R , mentioned earlier. In 
QCD (and quenched QCD), Goldstone bosons are 
emitted and absorbed by large, composite objects 
built of quarks and gluons. Whenever a compos- 
ite object emits or absorbs a probe with finite 
momentum one must have a form-factor which 
will suppress such processes for momenta greater 
than A ~ Rr 1 . Indeed, for > A we expect 
Goldstone boson loops to be suppressed as pow- 
ers of A/m„, not m T /A (or m„/ A x sb). Of course, 
this does not necessarily mean that one cannot in 
principle carry through the program of xpt. How- 
ever, it does mean that there may be considerable 
correlations between higher order coefficients and 
that it may be much more efficient to adopt an 
approach which exploits the physical insight we 
just explained. 

Over the past three years or so we have devel- 
oped an efficient technique to extrapolate every 
hadron property which can be calculated on the 
lattice from the large mass region to the physical 
quark mass - while preserving the most important 
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non-analytic behaviour of each of those observ- 
ables. This task is not trivial, in that various ob- 
servables need different phenomenological treat- 
ments. On the other hand, there is a unifying 
theme. That is, pion loops are rapidly suppressed 
for pion masses larger than A (ra^ > 0.4 — 0.5 
GeV). In this region the constituent quark model 
seems to represent the lattice data extremely well. 
However, for m w below 0.4-0.5 GeV the Gold- 
stone loops lead to rapid, non-analytic variation 
with vrtq and it is crucial to preserve the correct 
leading non-analytic (LNA) and sometimes the 
next-to-leading non-analytic (NLNA) behaviour 
of xpt. 

In order to guide the construction of an ef- 
fective, phenomenological extrapolation formula 
for each hadron property, we have found it ex- 
tremely valuable to study the behaviour in a par- 
ticular chiral quark model - the cloudy bag model 
(CBM) Built in the early 80's it com- 

bined a simple model for quark confinement (the 
MIT bag) with a perturbative treatment of the 
pion cloud necessary to ensure chiral symmetry. 
The consistency of the perturbative treatment 
was, not surprisingly in view of our earlier dis- 
cussion, a consequence of the suppression of high 
momenta by the finite size of the pion source (in 
this case the bag). Certainly the bag model, with 
its sharp, static surface, has its quantitative de- 
fects. Yet the model can be solved in closed form 
and all hadron properties studied carefully over 
the full range of masses needed in lattice QCD. 
Provided one works to the appropriate order the 
CBM preserves the exact LNA and NLNA be- 
haviour of QCD in the low mass region while nat- 
urally suppressing the Goldstone boson loops for 
TOtt > A ~ R . Finally, it actually yields quite a 
good description of lattice data in the large mass 
region. 

We now summarise the particular situation 
with respect to the chiral extrapolation of lat- 
tice data for some phenomenologically significant 
baryon properties. 

2. HADRON MASSES 

By far the most extensive and accurate data 
for hadron properties concerns their masses. At 



present the dynamical quark data is limited to 
quark masses a little over 5 times larger than the 
physical light quark masses, with most of the data 
at masses at least 10 times larger. One can expect 
to have data at 2-3 times the physical quark mass 
within 5 years and hence the crucial issues for 
chiral extrapolation are: 

• a) How accurately can one extract physical 
hadron masses from current data? 

• b) How accurately can one extract physical 
hadron masses from the next generation of 
data. 

• c) How model dependent are the results? 

Until the developments in Adelaide in the last 
few years p7[, the usual approach to chiral extrap- 
olation was to draw a straight line through the 
data. That is, taking the nucleon as the classic 
example: 



m N (m v ) = c + c 2 ml. 



(1) 



For pion masses above 1 GeV or so this does not 
give such a good fit. However, this can be at- 
tributed to the fact that rrvi is no longer propor- 
tional to m q in this region |14|,[l9|,^o| . We are not 
concerned with such large quark masses, instead 
we concentrate on the region m K < 1 GeV. Only 
recently has data in this region exhibited signifi- 
cant deviation from Eq. (Q) |2l],|22j and this led 
to the more sophisticated extrapolation function: 



m N {m v ) = c + c 2 ml + c 3 m^, 



(2) 



where the presence of the m\ term is motivated 
by chiral symmetry. (In terms of m q this term 

is non-analytic, being proportional m^ 2 , and is 
the LNA term in the formal expansion of the nu- 
cleon mass.) Indeed, it is relatively simple to see 
that C3 comes from the pion pole term in the 7rN 
self-energy process and hence that it is model in- 
dependent Q: 



C3 = clna 



3gj 
327r/a' 



(3) 



with gA and f n the axial charge of the nucleon 
and the pion decay constant, respectively, in the 
chiral limit. 
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The first empirical indication of serious prob- 
lems in this approach came with the realization 
that a fit to lattice data gives C3 ~ —0.76 GeV~ 2 , 
whereas the model independent value given by 
Xpt, Eq.(||), is -5.6 GeV -2 - a factor of 8 larger! 
This tells us immediately that either there are 
serious convergence problems with Eq. (g) or lat- 
tice QCD is in error. Clearly most readers would 
opt for the first possibility and so do we. 

The formal series for about the chiral 

SU(2) limit is usually written as 



mjv = mo + c 2 ml + clnaw4 + c ^ m \ 
+c NLN a"4. lnm* + c 5 ml + c 6 ml + .. 



(4) 



where the next-to-leading non-analytic term, 
m^lnm^, is dominated by the N — > Air —> N 
loop. However, for any meaningful extrapolation 
of lattice data in the next decade this expression 
is essentially useless. It is derived in the limit 
TOtt << A(= mi — m^), whereas the lowest lat- 
tice data with dynamical fermions that one can 
expect in the next decade is perhaps 200-250 MeV 
- c.f. A = 292 MeV. In fact, most lattice data 
will still lie above A. Mathematically the region 
around m n w A is dominated by a square root 
branch cut which starts at = A. Using di- 
mensional regularization this takes the form jlS| : 



25tt 2 / 2 
A 



(A 2 -m 2 )Mn(A 



(2A^ - 3m^)lnm w 



- \J A 2 - to 2 ) 
(5) 



for to w < A, while for to^ > A the first logarithm 
becomes an arctangent. No serious attempt has 
been made to extend the formal expansion in Eq. 
(Q) to incorporate this cut in an analysis of lattice 
data and, given the number of parameters to be 
determined if one works to order to^., it is not 
likely that it will be done in the next decade. 

Even ignoring the Air cut for a short 
time, studies of the formal expansion of the 
N — >N 7r — > N self-energy integral (ctntt, suggest 
that it has abysmal convergence properties. Us- 
ing a sharp, ultra-violet cut-off, Wright showed 
H that the series diverged for m w > 0.4 GeV. If 
one instead uses a dipole cut-off, which in view of 
the phenomenological shape of the nucleoli's axial 



form-factor is much more realistic, it is worse - 
with the radius of convergence being around 0.25 
GeV. 

In summary, the lattice data itself, formal stud- 
ies in the NJL model and phenomenological stud- 
ies of the nucleon based on long-distance regular- 
ization all suggest that the radius of convergence 
of the formal chiral expansion of the nucleon mass 
barely touches the lowest data point that we will 
get from lattice QCD in the next decade. This is 
not a productive approach to the problem. 

2.1. The Solution 

On the other hand, there is a better way. One 
can directly use the form [fl7|| : 

tojv = a + /3to 2 + o- N7r (m„., A) + cr A7T (m^,A), (6) 

where cb^ is the self-energy arising from a B it 
loop (B = N or A) and a, (3 and A are determined 
by fitting lattice data. For the reasons outlined, 
it is essential that the self-energies are evaluated 
using some ultra-violet regulator - a sharp cut- 
off or a dipole form, for example. Whatever is 
chosen does not effect the non-analytic structure 
which is guaranteed correct. The branch points 
at TOtt equals zero and A are incorporated natu- 
rally. While a and (3 will depend on the regula- 
tor, one can easily expand the self-energy terms 
to order to 2 (or higher) to obtain the chiral coef- 
ficients at the appropriate order to compare with 
effective field theory - e.g. see Ref. |l3| for a 
full discussion of this issue. The essential point 
is that studies of the nucleon, the A (c.f. Fig. 
4 of Leinweber et al. and the p meson |23|] 

suggest that this procedure will result in little or 
no model dependence in the extrapolation to the 
physical pion mass once there is accurate lattice 
data for to t ~ 0.3 GeV or less. Physically this is 
possible because the self-energy loops are rapidly 
suppressed in the region to t > 0.4 GeV. Thus, 
an extrapolation based on Eq. (^|) allows one to 
respect all the chiral constraints, keep the num- 
ber of fitting parameters low and yield essentially 
model independent results at the physical pion 
mass. No other approach can do this. 

2.2. A Possible Connection to QQCD 

The study of baryon spectroscopy in quenched 
lattice QCD (QQCD) has recently made great 
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progress. We have already noted that the lat- 
tice data behaves like a constituent quark model 
for quark masses above 50-60 MeV because Gold- 
stone boson loops are strongly suppressed in this 
region. This not only provides a very natural ex- 
planation of the similarity of quenched and full 
data in this region but it also suggests a much 
more ambitious approach to hadron spectra. It 
suggests that one might remove the small effects 
of Goldstone boson loops in QQCD (including the 
rf) and then estimate the hadron masses in full 
QCD by introducing the Goldstone loops which 
yield the LNA and NLNA behaviour in full QCD. 

As a first test of this idea, Young et al. [M re- 
cently analysed the MILC data |2^] for the N and 
A, using Eq.(^|) for full QCD and the appropriate 
generalization for QQCD - i.e. using quenched 
pion couplings as well as the single and double rf 
loops [ p5p^ |. The results were remarkable, with 
the values of a and (3 for the N (or the A) ob- 
tained in QQCD agreeing within statistical errors 
with those obtained in full QCD. Certainly this 
result is somewhat dependent on the shape of the 
ultra-violet cut-off chosen - although the extent 
of that is yet to be studied in detail. Nevertheless, 
given that the study involved the phenomenolog- 
ically favoured dipole form, it is a remarkable re- 
sult and merits further investigation. 

3. ELECTROMAGNETIC PROPERTIES 
OF HADRONS 

Although there is only limited lattice data for 
hadron charge radii, recent experimental progress 
in the determination of hyperon charge radii |^?| , 
has led us to examine the extrapolation procedure 
for extracting charge radii from the lattice simu- 
lations. Figure [l] shows the extrapolation of the 
lattice data for the charge radius of the proton, 
including the In (LNA) term in a generalised 
Pade approximant 
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<r 2 > p - 

^ ' ^ch 



C\ + X m 



1 + c 2 m£ 



(7) 



Here c\ and ci are parameters determined by 
fitting the lattice data in the large mass region 
(ra 2 > 0.4 GeV 2 ), while /j,, the scale at which the 



effects of pion loops are suppressed, is not yet de- 
termined by the data but is simply set to 0.5 GeV. 
The coefficient x is model independent and de- 
termined by chiral perturbation theory. Clearly 
the agreement with experiment is much better if, 
as shown, the logarithm required by chiral sym- 
metry is correctly included - rather than simply 
making a linear extrapolation in the quark mass 
(or to 2 ). Full details of the results for all the octet 
baryons may be found in Ref. pq] . 

The situation for baryon magnetic moments is 
also very interesting. The LNA contribution in 
this case arises from the diagram where the pho- 
ton couples to the pion loop. As this involves two 
pion propagators the expansion of the proton and 
neutron moments is: 



(8) 



fjPW=[/W T am n + 0(ml). 

Here /Uq is the value in the chiral limit and the 

i 

linear term in m w is proportional to rriq , a branch 
point at m q — 0. The coefficient of the LNA term 
is a = 4.4/XArGeV -1 . At the physical pion mass 
this LNA contribution is 0.6/ijv, which is almost 
a third of the neutron magnetic moment ]3C| ] . 
As for the charge radii, the chiral behaviour of 
vital for a correct extrapolation of lattice 
data. From Fig. |^ we see that one can obtain 
a very satisfactory fit to some rather old data 
(which happens to be the best available) using 
the simple Pade approximant 



M P(n) = 



p(n) 
Mb 



1 ± 



/3ml 



(9) 



Existing lattice data can only determine two pa- 
rameters and Eq.(^) has just two free param- 
eters while guaranteeing the correct LNA be- 
haviour (i.e. the correct value of a) as m T — ■> 
as well as the correct behaviour of HQET at 
large m 2 . The extrapolated values of fi p and 
/i™ at the physical pion mass, 2.85 ± 0.22/^jv and 
— 1.90 ± 0.15/iat are currently the best estimates 
from non-perturbative QCD (3l). (Similar re- 
sults, including NLNA terms in chiral perturba- 
tion theory, have been reported recently by Hem- 
mert and Weise p2|.) For the application of sim- 
ilar ideas to other members of the nucleon octet 
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Figure 1. Chiral extrapolation of lattice results 
for the squared electric charge radius of the pro- 
ton (solid line) - from Ref. pgfl . Fits to the 
contributions from individual quark flavors are 
also shown (the it-quark results are indicated by 
open triangles and the d-quark results by open 
squares). The physical value predicted by the fit 
is indicated at the physical pion mass, while the 
experimental value is denoted by an asterisk. 



Figure 2. Chiral extrapolations of lattice QCD 
magnetic moments for the proton (upper) and 
neutron (lower). The curves illustrate a two pa- 
rameter fit, Eq. (^|), to the simulation data, us- 
ing a Pade approximant, in which the one-loop 
corrected chiral coefficient of m ff is taken from 
Xpt. The experimentally measured moments are 
indicated by asterisks. The figure is taken from 
Ref. §l|. 



we refer to Ref. j33|, while for the strangeness 
magnetic moment of the nucleon we refer to Ref. 

. The last example is another case where there 
have been tremendous improvements in the ex- 
perimental capabilities. Specifically, the accurate 
measurement of parity violation in ep scattering 
is giving us vital information on hadron struc- 
ture. 

In concluding this section, we note that the ob- 
servation that chiral corrections are totally sup- 
pressed for m,r above about 0.4 GeV and that the 
lattice data looks very like a constituent quark 
picture there suggests a novel approach to mod- 
elling hadron structure. It seems that one might 
avoid many of the complications of the chiral 
quark models, as well as many of the obvious 
failures of constituent quark models by build- 
ing a new constituent quark model with u and 
d masses in the region of the strange quark - 
where SU(3) symmetry should be exact. Com- 



parison with data could then be made after the 
same sort of chiral extrapolation procedure that 
has been applied to the lattice data. Initial re- 
sults obtained by Cloet et al. for the octet baryon 
magnetic moments using this approach are very 
promising indeed |36| . We note also the extension 
to A-baryons, including the NLNA behaviour, re- 
ported for the first time in Ref. p7| . 

4. MOMENTS OF STRUCTURE FUNC- 
TIONS 

The moments of the parton distributions mea- 
sured in lepton-nucleon deep inelastic scattering 
are related, through the operator product expan- 
sion, to the forward nucleon matrix elements of 
certain local twist-2 operators which can be ac- 
cessed in lattice simulations [p| . The more recent 
data, used in the present analysis, are taken from 
the QCDSF |@ and MIT ||J groups and shown 



7 



a 0.04- 



QCDSF Cockeler et al., 1 
QCDSF Gockeler et al.. 1 
QCDSF Best et al., 1997 
MIT-DD60Q (Quenched) 
MlT-SliSAM (Full) 



Experimen 
Meson Cloi 
Model 




0.4 0.6 
m 2 [GeV 2 ] 



Figure 3. Moments of the u — d quark distribution 
from various lattice simulations. The straight 
(long-dashed) lines are linear fits to this data, 
while the curves have the correct LNA behaviour 
in the chiral limit - see the text for details. The 
small squares are the results of the meson cloud 
model and the dashed curve through them best 
fits using Eq. (|Tc|) . The stars represent the phe- 
nomenological values taken from NLO fits in the 
MS scheme. The figure is taken from Ref. 142], 



in Fig. |3| for the n = 1, 2 and 3 moments of the 
u — d difference at NLO in the MS scheme. 

To compare the lattice results with the exper- 
imentally measured moments, one must extrap- 
olate in quark mass from about 50 MeV to the 
physical value. Naively this is done by assuming 
that the moments depend linearly on the quark 
mass. However, as shown in Fig. || (long dashed 
lines) , a linear extrapolation of the world lattice 
data for the u—d moments typically overestimates 
the experimental values by 50%. This suggests 
that important physics is still being omitted from 



the lattice calculations and their extrapolations. 

Here, as for all other hadron properties, a lin- 
ear extrapolation in fa ~ m\ must fail as it omits 
crucial nonanalytic structure associated with chi- 
ral symmetry breaking. The leading nonana- 
lytic (LNA) term for the u and d distributions in 
the physical nucleon arises from the single pion 
loop dressing of the bare nucleon and has been 
shown [ftO 41 to behave as m^logm^. Expe- 
rience with the chiral behaviour of masses and 
magnetic moments shows that the LNA terms 
alone arc not sufficient to describe lattice data for 
?n > 0.2 GeV. Thus, in order to fit the lattice 
data at larger m^, while preserving the correct 
chiral behaviour of moments as — > 0, a low 
order, analytic expansion in rn\ is also included 
in the extrapolation and the moments of u — d are 
fitted with the form E^l: 



{x n ) u -d = a„ + b n ml 
+ a n clna ml In 



(10) 



where the coefficient, clna = — (3(7^ + l)/(4yT/ 7r ) 2 
f4l}| . The parameters a n , b n are determined by fit- 
ting the lattice data. The mass /j determines the 
scale above which pion loops no longer yield rapid 
variation and corresponds to the upper limit of 
the momentum integration if one applies a sharp 
cut-off in the pion loop integral. Consistent with 
our earlier discussion of this scale it is taken to 
be 0.55 GeV. Multi-meson loops and other con- 
tributions cannot give rise to LNA behaviour and 
thus, near the chiral limit, Eq. ( |To"| ) is the most 
general form for moments of the PDFs at 0(yn\) 
which is consistent with chiral symmetry. 

We stress that \i is not yet determined by the 
lattice data and it is indeed possible to consis- 
tently fit both the lattice data and the experi- 
mental values with /i ranging from 0.4 GeV to 
0.7 GeV. This dependence on /i is illustrated in 
Fig. |^ by the difference between the inner and 
outer envelopes on the fits. Data at smaller quark 
masses, ideally m\ ~ 0.05-0.10 GeV 2 , are there- 
fore crucial to constrain this parameter in order 
to perform an accurate extrapolation based solely 
on lattice data. 

Although we do not have sufficient space to 
explore the matter in detail, we note that once 
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we have a parametrization of the first four mo- 
ments as a function of pion mass it is possible to 
study the behaviour of the PDFs as a function of 
x and Tot,- |43| , Remarkably, for m„ > 0.5 GeV 
the distribution tends to peak at x ~ 1/3, resem- 
bling very much what one would expect in a con- 
stituent quark model. This further reinforces the 
discussion in sect. 1.2 concerning the behaviour of 
hadron structure within QCD as the light quark 
mass exceeds 50-60 MeV. 

4.1. Spin Dependent Moments 

There has been considerable interest in spin- 
dependent structure functions since the discovery 
of the spin crisis by EMC. As a result we now 
have a great deal of experimental information on 
spin-dependent PDFs . More recently studies 
of the axial charge of the nucleon within lattice 
QCD have revealed a surprisingly strong depen- 
dence on the volume of the lattice Q . Of course, 
interest is not limited to the first moment of the 
usual PDFs. One can calculate higher moments 
on the lattice, as well as moments of the trans- 
verse structure functions and experiments are also 
planned to explore the latter. 

The analysis of the chiral extrapolation of mo- 
ments of the unpolarized PDFs in the previous 
subsection has recently been extended to the po- 
larized sector by Detmold et al. jf6|. Years of 
phenomenological experience suggest that the A 
resonance might be important for spin-dependent 
observables. Within the CBM the explicit ap- 
pearance of the A played a vital role in the con- 
vergence of a perturbative expansion of nucleon 
properties, including its axial charge [Q. Thus 
it is not surprising that the explicit inclusion of 
the A in the extrapolation of the spin-dependent 
moments has a very significant effect. Whereas 
the rapid variation of the unpolarized moments 
is unaffected, the inclusion of the A almost to- 
tally removes the non-linearity for the polarized 
moments. While Detmold et al. suggest a new 
extrapolation formula which accurately approxi- 
mates the non-analytic behaviour in this case, one 
could fairly reliably use a linear chiral extrapola- 
tion for the spin-dependent moments. For a full 



discussion of the results we refer to Ref. 46 1, sim- 
ply noting here that with a careful extrapolation 



of data taken on a relatively large lattice 0] the 
agreement between theory and experiment for gA 
is now quite reasonable: ^(Latt) = 1.12 ± 0.05 
compared with (^(Expt.) = 1.267. 

5. CONCLUSION 

At the present time we have a wonderful con- 
junction of opportunities. Modern accelerator fa- 
cilities are providing data of unprecedented pre- 
cision over a tremendous kinematic range at the 
same time as numerical simulations of lattice 
QCD are delivering results of impressive accu- 
racy. It is therefore timely to ask how to use 
these advances to develop a new and deeper un- 
derstanding of hadron structure and dynamics. 

In combination with the very successful tech- 
niques for chiral extrapolation, which we have il- 
lustrated by just a few examples, lattice QCD will 
finally yield accurate data on the consequences of 
non-perturbative QCD. Furthermore, the physi- 
cal insights obtained from the study of hadron 
properties as a function of quark mass will guide 
the development of new quark models and hence 
a much more realistic picture of hadron structure. 
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